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Simulation of a Zenith-Field Sky Survey on RATAN-600 
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Abstract. In this paper we simulate a deep multi-frequency zenith-field sky survey on RATAN-600 (the RZF 
survey). In our simulations we use the 1.4-GHz sky images obtained in the NVSS survey. We convolved NVSS 
images with the two-dimensional power beam pattern of RATAN-600 and obtain simulated 24-hour scans of sky 
transits at all wavelengths of the RZF survey. For the 7.6-cm wavelength we analyze the effect of the image area 
size on the results of the simulation. We estimate the accuracy of the determination of source fluxes on simulated 
scans and derive the distributions of the spectral indices of the sources. We use the simulated scans to clean real 
records of the RZF survey at 7.6 cm. The standard error of the residual noise at this wavelength is about 1 mjy. 



1. INTRODUCTION 

Sky surveys are among the main sources of informa- 
tion about cosmic objects. The first deep sky survey on 
RATAN-600 was carried out in 1980-1982 during the 
"Cold" experiment at the declination of 5 ~ 5° pQ. The 
use of the best radiometer in terms of fluctuation sensi- 
tivity made it possible to obtain new constraints on the 
background radiation of the Universe and prepare a cata- 
log of radio sources with a detection threshold of ~ 10 mjy 
(RC catalog) at 7.6 cm [2j|3]. 

During the period from 1979 through 1987 the multi- 
frequency "Zelenchuk Survey" in the declination interval 
5 = 0° -7- 14° [1 |S], and during the period from 1987 
through 1989, a polar survey [5J, were carried out on the 
Southern sector with a periscopic reflector. The thresh- 
old of 3.9-GHz source detection flux was equal to 50 and 
14 mjy for the "Zelenchuk" and polar surveys, respec- 
tively. In 1988 the sky area 8 = 47°06'45" 47°07'45", 
8 h < a < 14 h was observed with the whole ring aperture 
of RATAN-600 with a detection threshold of ~ 15 mjy 
0. 

The next step in conducting surveys on RATAN-600 
was taken in the late 1990ies. The upgrade of the radio 
telescope and substantial improvement of the radiometer 
sensitivity made it possible to perform a deeper survey at 
a higher quality level. 

The RZF survey (RATAN-600 Zenith Field) was per- 
formed during the period from 1998 through 2003 on the 
Northern sector at wavelengths 1.0-^55 cm at declination 
of 384. At the main wavelength of observations (A7.6 cm) 
the half-power survey width was <520oo.o = 41°30'42" ± 2', 
h < CK2000.0 < As a result of the survey, a catalog of 
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radio sources at 7.6 cm (the RZF catalog) was compiled 
[8]. The minimum flux of the radio sources of the RZF 
catalog was close to about 2.5 mjy, which is comparable 
to the threshold flux of the NVSS catalog [£>]. Starting 
from 2001, eight sky bands were observed in addition to 
the central band. The additional bands were shifted by 
AS = ±12'n (n=l, 2, 3, 4) with respect to the central 
band. Note that the density of the survey bands in dec- 
lination was increased by a factor of two since 2006, and 
the interval between the neighboring bands became equal 
to 6'. 

A distinctive feature of the sky surveys performed with 
RATAN-600 compared, e.g., to VLA surveys (NVSS 0, 
FIRST [10]), is that at a given declination we observe 
sources located within a certain interval (band) of decli- 
nations. The angular size of this band is determined by 
the size of the vertical power beam (PB) of the radio tele- 
scope. The only exception is the sky survey performed in 
1988 in the mode of whole circular aperture of the antenna 
with a "pencil-shaped" beam of the radio telescope at the 
half-power level. In all other surveys one of the sectors of 
RATAN-600 was used (1/4 of the ring aperture). In this 
mode of operation the power beam pattern of the radio 
telescope differs substantially from that of a paraboloid. 

A distinctive feature of the power beam of RATAN-600 
operating in a single-sector mode is its large extension in 
the vertical plane and the extended "fan-shaped" back- 
ground at large elevations, whose intensity decreases as 
1/y in the vertical plane [HI [12] . Figure [T] gives an exam- 
ple of a PB computed for the elevation angle H — 88°42". 
Such a structure of the power beam allows, on the one 
hand, sources to be observed within a certain declination 
band and, on the other hand, makes the sources difficult 
to distinguish and thus results in extended structures ap- 
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Fig. 1. Power beam pattern of RATAN-600 computed for 
the wavelength and elevation of 7.6 cm and H = 88°42", 
respectively. 




Fig. 2. Fragment of the 4° x 4° NVSS image. The size in 
declination is trimmed (5 — ±20'). 

pearing on the records of sky transits across the fixed PB 
of the telescope (the so-called sky scans). This is because 
the magnitude and shape of the signal from an individual 
source depend substantially on the distance of the source 
from the central section of the power beam pattern. 

In view of the aforesaid, simulation of sky surveys con- 
ducted with RATAN-600 is of considerable practical in- 
terest. Such simulations are necessary both for studies of 
extended cosmic objects in order to estimate the contri- 
bution of point sources to their flux, and for extraction of 
numerous weak sources against noise background, thereby 
reducing the effect of "confusion" . 

In this paper we attempt to simulate the RZF sky sur- 
vey in its central section S a = 41°30'42" using the 1.4- 
GHz sky images (charts) obtained as a result of the NVSS 
(NRAO VLA Sky Survey) [9]. 

2. TECHNIQUE OF SIMULATION OF THE 
ZENITH-FIELD SURVEY 

Bursov and Majorova [13] were the first to simulate 
the zenith-field sky survey (RZF-survey). To this end, 



the above authors used the data of the NVSS catalog, 
from which they adopted the principal parameters of the 
sources. The program package for computing the two- 
dimensional PB of RATAN-600 [11] was used to simulate 
the curves of the transit of sources across various horizon- 
tal sections of the power beam pattern. These computa- 
tions assumed all sources to be points, i.e., their angular 
sizes were assumed to be much smaller than the halfwidth 
of the power beam (HPBW). In this case, the response of 
the antenna to the transit of a source with declination 8i S t 
across the power beam of the radio telescope is identical 
to the horizontal section of the power beam shifted by AS 
from the central section and multiplied by the 7.6-cm flux 
of the source in question. Here A<5 = Si St — <5o, where So is 
the declination of the central section of the RZF survey, 
which passed through 384. 

We computed the transit curves for most of the sources 
of the NVSS catalog, whose declinations fall within the 
AS = 41°30'42"±30' band. To compute the 7.6-cm fluxes 
of the sources, we used their 21-cm fluxes adopted from the 
NVSS catalog and the spectral indices of the correspond- 
ing sources. For the sources with unavailable spectral in- 
dices we set the spectral index equal to 7 = —0.8, which 
corresponds to a normal nonthermal spectrum (S u ~ v 1 ). 
We then time shifted the curves in accordance with the 
right ascension (a) of each particular source, and summed 
them up. The resulting one-dimensional scan of the diur- 
nal transit of the sky band across the fixed PB of the radio 
telescope provided a good analog of real observations. Let 
us refer this model of the RZF-survey as the model based 
on the NVSS catalog. 

In this paper we simulate the zenith-field survey by 
convolving the two-dimensional PB of RATAN-600 with 
the VLA NVSS sky images with a resolution of 9 = 45", 
which have the form of 4° x 4° areas. The images are in 
the FITS format and they are available from the NRAO 
site [TJ]. Each of the 2326 1024x1034 images has a name 
in the form PHHMMSDD, where P is the name of the 
Stokes parameter (I, Q, or U); and MM are the hours and 
minutes of right ascension of the image center; S, the sign 
of the declination (P and M for + and — , respectively), 
and DD is the declination of the image center in degrees 
for the epoch of 2000 (J2000). 

To simulate the RZF survey, we use about 20 m -long 
intensity images IhhOOPDD, Ihh20PDD, and Ihh40PDD 
(where hh = 00, 01, 23) with the central declination of 
S = +40° for the simulation at A < 7.6 cm; S = +40° 
and +44° for the simulation at A13 and 31 cm, and Sq = 
+36°, +40°, and +44°for the simulation at A49 cm. 

We compute the two-dimensional power beam pat- 
tern of RATAN-600 as described by Majorova pT|. We 
performed our computations for the elevation of 3C84 
(H = 88°42") in such declination intervals that maximum 
values of the power beam in the most distant sections 
are equal to a fraction of percent of the maximum of the 
diagram in its central section. The width of this declina- 
tion interval was equal to dS = ±6', dS = ±37'15", and 
dS = ±4°53' forXl.0, A7.6, and A49, respectively. We com- 
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Fig. 3. Normalized simulated 20-minute scans of the 62000.0 = 41°30'42", 12 h 30 m < a 2 ooo.o < 12 h 50 m sky area at 
wavelengths 1.0, 2.7, 7.6, 13, and 31 cm obtained by convolving the PB of RATAN-600 with the 4° x 4° areas of NVSS 
images. We normalized the scans to the signal level from 384. The curves are shifted only along the vertical axis. 



pute the power beam pattern of the radio telescope with a 
step equal to the pixel size on the NVSS survey maps (15 
arcsec). Figure Q] shows the two-dimensional power beam 
computed for the wavelength of 7.6 cm and elevation of 
H = 88°42". The size of the power beam in the figure is 
trimmed, dS = ±20'. 

We convolved the diagram with the part of the NVSS- 
survey images whose declination size is equal to the verti- 
cal size of the computed PB for the corresponding wave- 
length (dS). The center of NVSS images was located at 
the declination of 6 = 41°30'45", which virtually co- 
incides with the declination of 384. In addition to the 
IHHMMP40 images we used IHHMMP44 (at wavelengths 
7.6 cm ,Ai49 cm) and IHHMMP36 and IHHMMP44 (at the 
49-cm wavelength) images. Figure [5] shows an example of 
an NVSS image with a declination size of d5 — ±20'. 

We convert the images from the FITS format into 
binary matrices, which we then convolved with two- 
dimensional PB. As a result of these convolutions, we 
obtain 72 one-dimensional 20-minute sky scans at each 
wavelength. Note that the declination size of NVSS im- 
ages (and that of the power beam) exceed the width of the 
survey band at the wavelength considered by more than 
one order of magnitude, thereby allowing better account 
to be taken of strong distant sources falling within the PB 
of the radio telescope. Further increase of image sizes in 6 
has virtually no effect on the results of the simulation. 

We normalize the resulting scans to the level of the 
3C84 signal and convert them into F-format files. F format 
is a modification of FITS format used on RATAN-600 for 
representing the actual records and for their subsequent 
reduction. 

Figure [3] shows the normalized simulated 20-minute 
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12^50™ sky area at 1.0, 2.7, 7.6, 13, and 31 cm wavelengths 
obtained using the technique described above. Here and in 
subsequent figures scans are shifted with respect to each 
other along the vertical axis for better visualization. 

It is evident from the curves presented that as we 
pass from A1.0 cm to longer wavelengths the halfwidths 
of the sources in records increase and new sources and 
extended structures appear. This is due to the increase 
of the size of the PB with increasing wavelength, so that 
the power beam begins to cover the sky areas increas- 
ingly distant from the central section. Strong sources 
passing far from the central section in terms of 6 would 
produce extended features similar to the one we see in 
Fig. [3] at A13 and 31 cm in the right-ascension interval 
12 /l 42 m < a < 12 h A7 m . 

Figure [4] shows normalized simulated scans at 7.6, 13, 
and 31 cm (curves 1 ) and the corresponding normalized 
transit curves obtained within the framework of the RZF 
survey (curves 2). The real records shown in the figures 
were kindly provided by N. N. Bursov. They were obtained 
by averaging about 200 transits of the same sky area in 
the central band of the survey during the period from 
1998 through 2003 [IB] , The simulated and real scans are 
normalized to the level of the signal from 384. 

The subsequent analysis is based on the results of the 
simulation of the RZF survey at 7.6 cm. The sensitivity 
of the RZF survey is maximum at this wavelength and it 
is close to the sensitivity of the NVSS survey. 

It is of interest to compare the results of simulations 
not only with real records, but also with simulated scans 
based on the NVSS catalog [13] data. Figure [5] features 20- 
minute simulated sky scans obtained using NVSS images 
(curves 1 ) and simulated scans based on the data of the 
NVSS catalog (curves 3), as well as the corresponding 
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Fig. 4. Normalized simulated scans at the wavelengths of 
7.6, 13, and 31 cm obtained using NVSS images (curves 1 ), 
and normalized real averaged records of the same sky areas 
obtained within the framework of RZF survey (curves 2). 
The scans are normalized to the level of the 3C84 signal. 
The curves are shown shifted along the the vertical axis. 



averaged real records of the RZF survey (curves 2) at 
7.6 cm. All scans are normalized to the level of the 3C84 
signal. The simulated scans based on the data of NVSS 
catalog and averaged real records were kindly provided by 
N. N. Bursov. Records are shifted with respect to each 
other along the vertical axis. 

It is evident from a comparison of the curves shown 
in Fig. O that both models agree well with observational 
data, however, the small-scale structure of simulated scans 



obtained using NVSS images is closer to that of real sky 
records. The simulated scans based on the NVSS catalog 
are smoother as compared to curves 1 and 2. 

These differences are even more apparent in Fig. [6l 
which shows 10-minute portions of records. The differ- 
ences can be explained by the fact that simulated scans 
based on the data of the NVSS catalog contain only the 
NVSS sources in the given declination band that have 
fluxes above a certain level. At the same time, the sim- 
ulated scans constructed using NVSS images contain, in 
addition to NVSS sources, VLA detector noise convolved 
with the PB of RATAN-600, and the sources whose sig- 
nals are at the detection level and close to the noise of the 
NVSS survey. Such sources could have escaped the NVSS 
catalog. 

Note that the simulation of the RZF survey using the 
technique proposed by Bursov and Majorova [13] also has 
certain advantages. Thus the levels of the signals from the 
sources with known spectral indices are closer to the actual 
values for the model based on the NVSS catalog, because 
the fluxes of these sources were converted from A21 cm to 
A7.6 cm in accordance with their spectral indices. 

Such an individual correction of source fluxes is very 
problematic to perform in the model proposed in this 
paper. The fluxes can be converted only in bulk for all 
sources in the record using the average spectral index. 

3. SOME PROBLEMS THAT CAN BE SOLVED 
BY SIMULATING THE RZF SURVEY 

Consider now the problems that can be solved with the 
simulated scans of the RFZ sky survey obtained. 

Note that simulated scans based both on the data of 
the NVSS catalog and NVSS images have already been 
used to identify sources of the zenith-field sky survey and 
in the compilation of the RZF catalog [8]. The high de- 
gree of correlation between simulated and observed scans 
(~ 90%) allowed us to separate the contribution of dis- 
crete sources and background radiation, and to identify 
the sources in the survey records with the sources of the 
NVSS catalog. Simulated curves are especially valuable 
for detecting weak objects and for finding new objects un- 
detected at decimeter-wave frequencies. 

In this paper we also address a number of other prob- 
lems. First, the analysis of simulated scans based on NVSS 
images allows us to estimate the accuracy of identification 
of sources and of the determination of their fluxes on sky- 
transit records. Second, simulated scans can be used to 
clean a real record from sources and estimate the mini- 
mum residual noise. Furthermore, it is of interest to esti- 
mate the spectral indices of the sources from the ratio of 
their signals on real and simulated scans and to compare 
their distribution with the distribution of spectral indices 
computed by Bursov et al. [5]. 

Before we start addressing these problems, let us list 
the main factors that may affect the accuracy of the sim- 
ulation based on NVSS images. These factors include: the 
accuracy to which we know the power beam pattern of the 
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Fig. 6. Ten-minute portions of normalized simulated sky scans obtained using 4° x 4° NVSS images (curves 1); the 
corresponding averaged records of the RZF survey (curves 2) at 7.6 cm, and simulated scans based on the data from 
the NVSS catalog [13]. 



radio telescope; the correctness of the source extraction on 
the simulated scans obtained, and the effect of the sizes of 
NVSS image areas on the results of simulation. 

The values of the vertical PB are used to convert the 
antenna temperatures of the sources to the corresponding 
antenna temperatures for the central section of the survey 
and for the subsequent determination of the source fluxes. 
Majorova and Trushkin [T^] and Majorova and Bursov 
[T7] show that if the antenna is in good state, then the 
computed PB agrees with the experimental PB very ac- 
curately and therefore the errors due to the use of the 
computed PB are minimal and amount to mere 1 -f- 3%. 

Extraction of sources on simulated and real scans us- 
ing the technique of Gauss analysis involves certain dif- 
ficult aspects, especially when it comes to weak sources 
and sources with close right ascensions. We estimate the 
accuracy of the determination of the fluxes of the sources 
identified on records in Paragraph 6. 

As for the effect of the size of the NVSS images on 
the results of simulation, our main task is to determine 
the maximum size of the area to minimize the simulation 
errors. 



4. EFFECT OF THE SIZE OF NVSS IMAGES 
ON THE SIMULATION RESULTS 

To perform the task mentioned above using simulated 
scans based on 7.6-cm NVSS images, we identify the 
sources of the NVSS catalog located within the ±6'-band 
of the central declination of the survey. Like in the case of 
real records, we use the Gauss analysis [16] to identify 
the sources, and control the declinations of the sources by 
the halfwidths of the Gaussians. To this end, we use the 
computed dependences of the halfwidth of the PB on its 
distance from the central section. 

We identified a total of about 500 NVSS catalog 
sources with 21-cm fluxes S > 2.7 mJy. We determined 
the signal level for each of the sources identified both in 
relative units (in terms of the 3C84 signal level) and ab- 
solute units (mJy). 

We construct the AS and Aa dependences of S m (0)/Sk 
for this sample of sources and determine the vertical PB of 
RATAN-600 and the dH dependence of the halfwidth of 
the power beam. Here Sk is the flux of the source according 
to the NVSS catalog; 5 m (0), the flux of the same source 
as inferred from the simulated scan and reduced to the 
central section of the survey (AS = 0), and dH is the 
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Fig. 7. Vertical PB constructed from sources extracted on simulated scans (filled squares) for the entire sample of 
sources (S > 2.7 mJy) (a) and from sources with fluxes S > 200 mJy — (b). Simulated scans were obtained using 
4° x 4° NVSS images. The solid lines indicate the vertical section of the computed PB at the wavelength of 7.6 cm 
after its convolving in each horizontal section with a Gaussian of width 0q.5 = 45". 
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Fig. 8. The dH dependences of the halfwidth of PB constructed using the data for the sources extracted on simulated 
scans (filled squares) for the entire sample of sources (S > 2.7 mJy) — (a) and using sources with fluxes S > 200 mJy — 
(b). Simulated scans were obtained from 4° x 4° NVSS images. The solid lines show the halfwidths of the computed 
PB at the wavelength of 7.6 cm convolved in each horizontal section with a Gaussian with a width of 9q.5 — 45". 
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Fig. 9. The AS dependences of S m (0)/Sk constructed using the entire sample of sources extracted on simulated scans 
(S > 2.7 mJy) (a) and using sources with fluxes S > 30 mJy (b) and S > 200 mJy (c). The simulated scans were 
obtained from 4° x 4° NVSS images. 



elevation offset with respect to the central section of the 
PB. 

S m (0) = S m (AS)/F r (dH), (1) 

where S m (AS) is the flux of the source with the coor- 
dinates (Sist, oust) as inferred from the simulated scan; 
F r (dH), the vertical PB computed using computer codes 
in accordance with the technique described by Majorova 
dJ, AS = S ist - S , Aa = a ist - a . 



During the survey considered the central section of the 
PB passes through the central section of the survey, im- 
plying that dH = AS. 

The values F of the vertical PB of the radio telescope 
for various dH can be inferred from the following formula: 

F(dH) = S m (AS)/S k , (2) 

where dH = AS, using the source fluxes inferred from 
simulated scans and the fluxes of the same sources listed 
in the NVSS catalog. 
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Fig. 10. The Aa dependences of S m (0)/Sk constructed using the entire sample of sources identified on simulated 
scans (S > 2.7 mJy) (a) and using sources with fluxes S > 30 mJy (b) and S > 200 mJy (c). The simulated scans 
were obtained from 4° x 4° NVSS images. 
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Fig. 11. Normalized five-minute simulated sky scans obtained using 1° x 1° NVSS images (curves 1) and the corre- 
sponding averaged real records of the RFZ survey (curves 2) at the wavelength of 7.6 cm. All scans are normalized to 
the 384 signal level. The curves are shown shifted along the vertical axis. 



The halfwidths of the power beams are equal to the 
halfwidths of the Gaussians fitted when identifying sources 
on simulated scans. 

The filled squares in Figs. [71(a) and[U(a) show the ver- 
tical power beams and the dH dependences of HPBW 
constructed from the entire sample of sources identified 
on simulated scans, and Figs. \T[b) and [Sfb) show the 
same power beams and dependences, but constructed from 
bright sources with fluxes S > 200 mJy. The solid lines 
show the vertical PB and the computed dH dependence of 
HPBW. Note that the two-dimensional power beams were 



convolved with the images already convolved with the in- 
strumental function of the VLA radio telescope, and there- 
fore we convolved the computed PB of RATAN-600 in each 
of the horizontal sections with a Gaussian of halfwidth 
equal to O .5 = 45" (HPBW VLA 0). 

A sufficiently strong discrepancy between the power 
beams constructed from simulated scans and the com- 
puted curves is immediately apparent from Figs. [7j and 
|U The pattern of the discrepancy does not depend on 
the brightness of the sources used to construct the depen- 
dences in question. 
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Figure [H shows the AS dependences of S m (0)/Sk con- 
structed for the entire sample of sources 
(S > 2.7 mJy) (a) and for the subsamples of sources with 
fluxes S > 30 mJy (b) and S > 200 mJy (c). All these 
dependences have a well-defined slope, indicating that the 
source fluxes inferred from simulated scans differ from the 
fluxes of the same sources listed in the NVSS catalog, and 
the discrepancy increases with the difference between the 
declinations of the source and of the central section of 
the survey. And, finally, the Aa dependences of S m (Q)/Sk 
shown in Fig. [TO] demonstrate that the discrepancy be- 
tween S m (0) and Sk increases with increasing distance 
between the source and the scan center (or between the 
source and the centers of image areas Ihh20P40, Ihh00P40, 
and Ihh40P40) . These features are especially apparent for 
the dependences based on the data for strong sources with 
fluxes S > 30 mJy and S > 200 mJy. 

All these discrepancies are due to projection of the sky 
points (i.e. points of a sphere) onto the planes of the NVSS 
images. Condon et al. [9] give the formulas that relate 
celestial (a, S) coordinates to pixel coordinates (x,y) of 
the image area: 

x = Xq — cosSsin(a — ao)/e, (3) 

y = Ho + [sinSsinSo — cosSsinSocos(a — ao]/e, (4) 

where e = tt/43200 radian = 15", x Q = 512, y = 513. 

It is evident from the above formulas that the farther 
is the source from the center of the sky area (ao, <5o)j the 
greater is the error of its inferred position in the image 
plane. The difference between values of AS and Ay = 
Hist ~ yo results in systematic errors in the inferred source 
fluxes reduced to the central section of the survey using 
the vertical power beam of the telescope. These errors, in 
turn, result in non-zero slopes of the dependences shown 
in Fig. [9] deviation of the simulated data points from the 
computed curves in Figs. [7] and and in the increase of 
the scatter of data points toward the scan boundaries in 
Fig. [TO] 

The dependences shown in Fig.[TOJcan be used to select 
the maximum right-ascension sizes of image areas such 
that the S m (0)/Sk ~ 1 condition is satisfied throughout 
the entire length of the simulated scan. This size is on the 
order of da = ±2.5™ (or ~ 1°). The declination size of 
the images, which is used in simulations, is close to its 
optimum value. Thus for the wavelength of A7.6 cm it is 
dS = ±37'15". 

Thus the simulated scans obtained by convolving 20- 
minute (in right ascension) sky areas of the NVSS survey 
with the two-dimensional PB of RATAN-600 are of great 
interest for the extraction of sources on real records of 
the zenith-field survey. The simulated and real scans are 
highly correlated and this correlation is especially useful 
for identifying weak sources. However, they are not suit- 
able for quantitative estimates. To this end, the size of the 
image areas must be on the order of 1° x 1° (about 5™ in 
terms of right ascension). Note that with such an image 
size both the number of images for the 7.6-cm wavelength 



and the number of simulated scans increases by a factor 
of four. This was one of the reasons why we began our 
simulations with 4° x 4° NVSS image areas. 



5. SIMULATION OF THE ZENITH-FIELD 
SURVEY USING 1° x 1° NVSS IMAGES 

We simulated the zenith-field survey using 1° x 1° NVSS 
images at the wavelength of 7.6 cm. We downloaded the 
images of the sky areas in FITS format from the site of 
SkyView [15] and converted them into the F format and 
then into 300x300 binary matrices. We then convolved 
the latter with the two-dimensional PB of RATAN-600 
computed for A7.6 cm and having the same mesh size in 
a and S as the images of the NVSS survey. 

The centers of the image areas had the following co- 
ordinates: a a = H h 5 m n, S = +41°30'45" + A, where H 
are hours from to 23, n — 0,1, 2, ...11, and A is the 
precession correction. 

As a result, we obtained a total of 288 five-minute sim- 
ulated scans. Figure. [TT1 (curves 1) shows examples of such 
scans. The same figure shows the averaged five-minute 
records of sky transits obtained in the RFZ survey at the 
wavelength of 7.6 cm during the 1998-1999 period (curves 
2). The simulated and real scans shown in the figures are 
normalized to the 384 signal level. The curves are shifted 
with respect to each other along the vertical axis. 

As is evident from the figures, the simulated and ob- 
served scans are highly correlated. The correlation coeffi- 
cient is ~ 0.95, which is higher than the correlation coeffi- 
cient between 20-minute simulated scans and the scans of 
the RFZ survey, and than the correlation coefficient be- 
tween the simulated scans based on the data of the NVSS 
catalog and real records. 

We use Gaussian analysis to identify the sources of 
the NVSS catalog on simulated scans; determine the am- 
plitude (in relative and absolute units) and halfwidth of 
their signals. We then use the results obtained to construct 
the same dependences as those constructed for the sources 
identified on 20-minutc simulated scans, namely: the AS 
and Aa dependences of S m (0)/Sk, the dH dependence of 
HPBW, and the vertical PB of the radio telescope. We 
show the results in Figs. [TOJ [TO] [TH and [TO] 

As is evident from the figures, the vertical PB and the 
dH dependence of HPBW constructed using five-minute 
simulated scans agrees well with the computed curves; the 
AS dependences of S m (0)/Sk have zero slopes, and the Aa 
dependences of S m (0)/Sk show a uniform scatter of data 
points along the entire scan length. This fact indicates that 
simulated scans contain no systematic errors and that they 
can be used to obtain qualitative estimates and not just 
to detect and identify the sources of the NVSS catalog on 
real scans. 
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Fig. 12. Vertical PB constructed using the sources extracted on simulated scans (filled squares): for the entire sample 
of sources (S > 2.7 mJy) (a) and for sources with fluxes S > 200 mJy (b). The simulated scans were obtained from 
1° x 1° NVSS images. The solid lines show the vertical section of the computed PB at 7.6 cm after its convolving in 
each horizontal section with a Gaussian of size 0q.5 — 45 . 
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Fig. 13. The dH dependences of the halfwidth of the vertical PB constructed using the sources identified on simulated 
scans (filled squares): for the entire sample of sources (S > 2.7 mJy) (a) and for sources with fluxes S > 200 mJy 
(b). The simulated scans were obtained from 1° x 1° NVSS images. The solid lines show the vertical section of the 
computed PB at 7.6 cm after its convolving in each horizontal section with a Gaussian of size #0.5 = 45 . 
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Fig. 14. The AS dependences of S m (0)/Sk constructed using the entire sample of sources identified on simulated 
scans (S > 2.7 mJy) (a) and using sources with fluxes S > 30 mJy (b) and S > 200 mJy (c). The simulated scans 
were obtained from 1° x 1° NVSS images. 



6. ESTIMATION OF THE ACCURACY OF THE 
INFERRED SOURCE FLUXES 



Identification of NVSS sources on simulated scans ob- 
tained by convolving NVSS images with the power beam 
pattern of RATAN-600 and subsequent reduction of their 
fluxes to the central section of the survey must yield the 
fluxes given in the NVSS catalog. However, the depen- 
dences shown in Figs. ITHfa.b.c) and [15](a,b,c) demon- 



strate evident scatter of data points about S m (0)/Sk = 
1.0, and this scatter increases with decreasing source flux. 

This is because as source fluxes decrease, source sig- 
nals and the noise from background sources become close 
in magnitude, thereby making extraction of sources diffi- 
cult both on simulated and real scans. Cases occur where 
sources blend together (are superimposed on each other). 
In some cases even strong sources are identified with insuf- 
ficient confidence if they have close right ascensions, but 
are located at different distances from the central section 
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Fig. 15. The Aa dependences of S m (0)/Sk constructed using the entire sample of sources identified on simulated 
scans (S* > 2.7 mJy) (a) and using sources with fluxes S > 30 mJy (b) and S > 200 mJy (c). The simulated scans 
were obtained from 1° x 1° NVSS images. 




Fig. 16. Histograms of spectral indices obtained using the fluxes listed in the NVSS catalog, 7 = Iu(St.q/ S21) (a) 
(open bars), and using the fluxes inferred from simulated scans based on NVSS images, 7 = ln(Sr.e/ $21) 3 ) (open 
bars). The histogram of spectral indices based on the data of the RZF catalog (filled bars) [5]. The histograms are 
based on a total sample of 317 sources with allowance for errors of inferred spectral indices. 



of the survey. All the above effects result in the distortion 
of signals from sources and even in the appearance of false 
sources. 

The errors of the inferred fluxes can be estimated using 
the resulting Aa dependences of 5 m (0)/5fc. The standard 
error cr™ is equal to the standard deviation of S m (Q)/Sk 
from its mean value averaged over the entire sample of the 
sources identified. 

This error gives the lower limit for the error of the de- 
termination of fluxes of sources identified on real scans at 
A7.6 cm. Recall that the error of the inferred fluxes on real 
scans also includes the error of interpolation (extrapola- 
tion) of fluxes to the wavelength of 7.6 cm. Other error 
components should be of similar magnitude, because we 
use the same method to identify sources and use the same 
instrumental function of the radio telescope to convert the 
source fluxes to their values in the central section of the 
survey. 

The standard error of the determination of the fluxes 
of the sources identified on simulated scans is: a™ = 
0.20 ± 0.01 — for the entire sample of the sources con- 
sidered with fluxes S > 2.7 mJy; cr™ = 0.10 ± 
0.01, for sources with S > 30 mJy and erg 1 = 0.08 ± 0.01, 
for sources with S > 200 mJy. (Here we give the source 
fluxes at A21 cm). 



Let us now compare the errors of cr<f with the errors 
a r s of the determination of the fluxes of sources identified 
in real records. The errors a r s determined by Bursov et 
al. [8] and Majorova and Bursov [17] based on the entire 
sample of sources considered exceed ex™ by 20-^25% (<r r s = 
0.24±0.03 pj], 0s = 0.25 [8]). For strong sources the error 
a r s was equal to 0.15 ± 0.03 [51Q7]. 

As we pointed out above, such discrepancies between 
<7™ and a r s are due primarily to errors of interpolation (ex- 
trapolation) of fluxes from A21 cm to A7.6 cm. However, 
a certain deviation of the real PB from the computed PB 
during the survey may also contribute to the above dis- 
crepancy. 



7. DETERMINATION OF THE SPECTRAL 
INDICES OF THE SOURCES 

The spectral index of a source can be computed as the 
logarithm of the Sr.e/S^i ratio or the logarithm of the 
Srs/S^l ratio, where S7.6 is the flux of the source as in- 
ferred from the real averaged record of the RZF survey at 
A7.6 cm; is the A21 cm flux of this source according to 
the data of the NVSS catalog, and S%{ is the source flux 
in mJy as inferred from the simulated scan. 
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Given that we use the same technique to identify 
sources on real and simulated scans, and given the high 
degree of correlation between these scans, one may as- 
sume that the errors of the extraction of the same source 
should be the same for both scans. That is, we should 
either overestimate or underestimate the flux of a source 
when identifying it simultaneously in the simulated and 
real records. Thus the spectral indices computed by the 
formula j m = ln(Si.§/ S^i™) may prove to be even more 
accurate than those inferred from the data of the NVSS 
catalog (jk — ln(Sr .6 / 'S21)) ■ This is true if the computed 
and real power beams of the radio telescope agree well 
with each other. 

To compare the spectral indices of the sources obtained 
using different methods, we identify on the averaged 7.6- 
cm records of the RZF survey the sources of the NVSS cat- 
alog within the ±6' band of the central section of the RZF 
survey. We identified a total of about 500 NVSS-catalog 
sources on the 24-hour record of the sky scan, determine 
the spectral indices jk and 7 m for each of these sources, 
and compare them to the spectral indices obtained by 
Bursov et al. [8]. For a more correct comparison, we used 
only the sample of 317 sources in common with the RFZ 
catalog. 

The open bars in Fig. \W\ show the histograms of the 
distribution of spectral indices obtained using the fluxes 
listed in the NVSS catalog, 7fc = ln(St ,% / (a) and the 
distribution of spectral indices computed using the fluxes 
obtained from simulated scans 7 m = ln(Sr.e/ S21) (b). The 
histograms are based on the sample of 317 sources that 
we identified on simulated and real scans. The filled bars 
in the same figure show the histograms of the spectral 
indices of the same sample of sources based on the data of 
the RZF catalog [8]. These histo grams were constructed 
with allowance made for the errors of the determination 
of spectral indices. 

As is evident from the figures shown, the histogram of 
the spectral indices j m computed using the fluxes inferred 
from simulated scans proves to be closer to the distribu- 
tion of spectral indices according to the data of the RZF 
catalog for the source sample considered. Given the con- 
fidence interval, they can be considered to be virtually 
identical. 

It follows from this that given the real records and 
the corresponding simulated scans obtained using NVSS 
images, one can infer the spectral indices of the sources 
from the ratio of their signals in the scans considered. The 
real and simulated scans must be normalized prior to that. 
Such a procedure is very easy and convenient to perform 
for sources whose spectral indices are unknown. 

8. ESTIMATE OF RESIDUAL NOISE ON REAL 
RECORDS USING SIMULATED SCANS 

We use the simulated scans based on NVSS images to 
clean the real records from discrete sources. To this end, 
we "join" the normalized five-minute scans into a single 



24-hour scan and convert it into a 7.6-cm wavelength scan 
in accordance with the average spectral index. 

As we pointed out above, in this method it is hardly 
possible to correct the signal from each source in accor- 
dance with its own spectral index even if the spectral in- 
dex is known, and therefore we used for our conversion the 
spectral indices close to those of the sources with normal 
nonthermal spectra, which lie in the 7 = 0.78 ± 0.02 in- 
terval. We chose the spectral index so as to minimize the 
standard error of the deviation of the simulated record cor- 
rected for the spectral index from the real record. We used 
trimmed records h < a < 2Q h (with the Galactic plane 
excluded) to estimate the magnitude of residual noise. 

The standard error of residual noise at 7.6 cm was 
equal to a — 7.6±0.9 mJy. This residual noise includes the 
noise due to unresolved sources; radiometer noise; atmo- 
spheric noise, and residual signals from the sources whose 
spectral indices differ from the average spectral index used 
to convert the source fluxes. 

Subtraction of the background inferred with a 30 s 
smoothing window (hereafter referred to as the thirty- 
second background) partially eliminates atmospheric noise 
and the noise from the sources located outside the AS = 
41°30'42" ± 6' band. This procedure reduces a of residual 
noise down to 4.6 ± 0.6 mJy. Note, for comparison, that 
the magnitude of residual noise after we subtract from the 
real averaged record the simulated scan obtained from the 
data of the NVSS catalog and then the thirty-second noise, 
is equal to a = 5.9 ± 0.6 mJy. This means that cleaning 
real records using simulated scans obtained by applying 
two different methods yields similar results. 

Figure IT7T a) shows the residual noise on a real aver- 
aged record of the RZF survey after subtracting from it the 
simulated scan based on NVSS images and corrected for 
the average spectral index. Figure [T7T b) shows the resid- 
ual noise obtained after further subtraction of the back- 
ground computed with a 30 s smoothing window from the 
difference of the real and simulated scans. 

As is evident from the figures, the records still contain 
signals from the sources whose spectral indices differ from 
the average spectral index. After eliminating such signals 
(they can be easily seen in the records) the standard devia- 
tion a of the residual noise becomes equal to about 1 mJy. 
Figure [T7T c) shows the record of the residual noise. Note 
that the standard error of the same portion of the un- 
cleaned real record (0 /l < a < 20 h ) is equal to 19 ±3 mJy. 
We did not take the 3C84 signal into account when esti- 
mating a. Thus cleaning of records using simulated scans 
reduces the residual noise of a real record by more than 
one order of magnitude. 

Note that records can be cleaned even without the use 
of simulated scans, simply by "cutting" the sources on 
records in accordance with their coordinates. However, 
such a cleaning procedure is very time consuming and 
sometimes ambiguous, because the transit curves obtained 
in RATAN-600 surveys are superpositions of source sig- 
nals convolved with the PB at different declinations. The 
proposed method allows easy elimination of most of the 
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point sources within the given declination band and with 
fluxes greater than 3 mJy (> 3<r). More refined methods 
are required for a more in-depth cleaning, and we do not 
describe them here, because they are a subject of a sepa- 
rate paper. 

9. CONCLUSIONS 

We propose and develop a method for the simulation of a 
deep multi-frequency zenith-field sky survey by convolv- 
ing NVSS images with two-dimensional power beams of 
RATAN-600. 

We perform simulations at 1.0, 2.7, 3.9, 7.6, 13, 31, and 
49 cm in the central band of the survey. We obtained 24- 
hour simulated scans of the sky transit across the power 
beam pattern of the telescope at each of the wavelengths 
listed above and compared the simulated scans with real 
records and with the simulated scans based on the data 
of the NVSS catalog [13]. 

We analyze the effect of the size of NVSS-image areas 
on the results of simulations at 7.6 cm. We show that simu- 
lated scans based on 4° x4° (20 m in right ascension) NVSS 
image areas are highly correlated (with a correlation coef- 
ficient of ~ 90%) with real scans and can therefore be used 
to identify and extract weak sources in the records of the 
RZF survey. However, for quantitative estimates the size 
of NVSS image areas should be reduced down to 1° x 1° 
(5 m in right ascension). 

We obtained a total of 288 five-minute simulated scans 
of sky transits at the wavelength of 7.6 cm by convolving 
1° x 1° NVSS images with the PB of RATAN-600, and then 
"joined" them into a single 24-hour scan. We identified a 
total of about 500 NVSS-catalog sources simultaneously 
in simulated scans and real averaged records of the RZF 
survey within the ±6' band of the central section of the 
survey and estimated their fluxes. 

The data obtained allowed us to estimate the accu- 
racy of the determination of source fluxes on simulated 
scans. The standard error of the determination of source 
fluxes on simulated scans is = 0.20 ± 0.01 for the entire 
sample of sources considered (S > 2.7 mJy). This error is 
20 -j- 25% lower than the standard error for real records 
0[T7]. The standard error is a™ = 0.10 ± 0.01 for sources 
with fluxes exceeding 30 mJy The standard errors <r™ ob- 
tained can be viewed as a lower limit for the errors of the 
determination of source fluxes on real scans at A7.6 cm. 

We derived the distributions of spectral indices of the 
sources by identifying sources on real and simulated scans. 
To compute the spectral indices, we used the source fluxes 
listed in the NVSS catalog and the source fluxes inferred 
from real and simulated scans. We compared the distribu- 
tions of spectral indices obtained using different methods 
with the distribution of spectral indices based on the data 
of the RZF catalog for the common sample of 317 sources. 
We show that the distribution of spectral indices com- 
puted as the logarithm of the ratio of the fluxes of sources 
identified on real and simulated scans agrees within the 
quoted errors with the distribution of spectral indices in- 



ferred from the data of the RZF survey of the same sample 
of sources. 

We use the simulated scans obtained to clean the real 
20-hour record of the RZF survey. After the subtraction of 
the thirty-second noise and elimination of residual signals 
of the sources whose spectral indices differ from the mean 
spectral indices the standard deviation a of residual noise 
was equal to ~ 1 mJy at 7.6 cm. This procedure reduced 
the noise of the initial record of the RZF survey by more 
than one order of magnitude. 
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Fig. 5. Normalized 20-minute simulated sky scans ob- 
tained using 4° x 4° NVSS images (curves 1); the cor- 
responding real records of the RZF survey (curves 2) at 
A7.6 , and simulated scans based on the data of NVSS cat- 
alog [13] (curves 3). All scans are normalized to the 3C84 
signal level. The curves are shown shifted with respect to 
each other along the vertical axis. 
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Fig. 17. Residual noise on a real averaged record of the 
RFZ survey after subtracting the simulated scan obtained 
using NVSS images (a); residual noise after subtracting 
the background obtained with a 30 s smoothing window 
from the difference of the real and simulated scans (b) , and 
residual noise after subtracting the sources whose spectral 
indices differ from the average spectral index () . 



